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Abstract. A shifted or misaligned feedback loop gives rise to a two-point nonlocality that is the spatial 
analog of a temporal delay. Important consequences of this nonlocal coupling have been found both in 
diffusive and in diffractive systems, and include convective instabilities, independent tuning of phase and 
group velocities, as well as amplification, chirping and even splitting of localized perturbations. Analytical 
predictions about these nonlocal systems as well as their spatio-temporal dynamics are discussed in one 
and two transverse dimensions and in presence of noise. 

PACS. 42.65 Sf Dynamics of nonlinear optical systems; optical instabilities, optical chaos and complexity, 
and optical spatio-temporal dynamics - 89.75.Kd Patterns 42.55.-fLasers 

1 Introduction domain if it is not reflected at the boundary. The exis- 

tence of these instabilities has been shown in systems in 



There is a considerable interest in dynamical regimes in 
which small fluctuations and "noise" are amplified. In a 
large class of optical systems this behavior is caused by 
convective instabilities |1|2| . A convective instability hap- 
pens in presence of some source of drift when a state of a 
nonlinear system becomes unstable and the group veloc- 
ity of a localized perturbation is larger than the velocity It was recently shown that significantly larger windows 



which the propagation of disturbances is characterized by 
drift or walk-off, modelled by a gradient term. For such 
systems small regions of convective instabilities have been 
predicted and observed in hydrodynamics 131, plasma [1] 
physics, traffic flow 5 and optics |1I6I7| . 



of propagation of the instability front. As a result, in the of convective instabilities are induced by nonlocal terms in 
laboratory frame the perturbation is amplified but moves the governing equations [8^. In optics, these terms result 
away and eventually decays at any point in a finite spatial from the presence of an off-axis or shifted feedback loop 



R. Zambrini, F. PapofF: Nonlocal feedback in nonlinear systems 




Fig. 1. Schematic representation of a nonlinear device with 
shifted feedback. The output beam is re-injected into the device 
after a shift Ax. 

which is modeled by a two point nonlocality that is the 
spatial analogous of a temporally delayed feedback. This 
is a common experimental issue and has been subject of 
both theoretical and experimental study in liquid crys- 
tals light valves |9ll0lllj . Kerr-like media [1216113) and 
generic nonlinear systems with diffusive [8^ and diffractive 
[T^ coupling. We note that the importance of feedback 
loops goes well beyond the realm of optics [15] and has 
been long recognized in other fields of physics and also 
in biology and engineering [16117] . Usually feedback loops 
are introduced to better control the system and to limit 
the growth of noise, while in the papers quoted a nonlo- 
cal feedback has been used mainly to study fundamental 
properties of fluctuations in non linear systems when con- 
vective instabilities are induced. In a recent experiment, 
however, an off-axis feedback loop has been used very ef- 
fectively to suppress noise-sustained structures caused by 
an intrinsic drift term in a free electron laser [TB], show- 
ing that two-point nonlocality are not only fundamentally 
interesting, but also extremely useful. 



In this paper we review the effect of off-axis feedback 
in a broad class of optical devices and explore the effect of 
the two-point nonlocality in the case of two transverse di- 
mensions. We consider passive as well as active nonlinear 
media with fast decay of the polarization [19| . including 
media with negative refraction index [20121] and devices 
with soft apertures [22]. Within the media we consider 
both diffusion and diffraction and show how an off-axis 
feedback changes the nature of the first instability thresh- 
old. As a consequence, there are large windows of control 
parameters where small localized signals can be strongly 
amplified while the background radiation in other region 
of the system remains very low [14] . The amplification 
does not need a continuous signal injection and takes place 
even when the initial perturbation is switched off. Further- 
more, in systems with diffraction and active media, the 
signal moves across the cavity with transverse phase and 
group velocities that are easily managed to have the same 
or opposite signs [H] . In spite of the broken transverse re- 
flection symmetry, localized perturbations can move both 
towards or against the off-set direction and can even split 
into two counter-propagating components, with the laser 
operating as a signal splitter. Both noise sustained struc- 
tures and signals control are shown by numerical simula- 
tions of the full nonlinear model confirming our theoret- 
ical analysis and a rich spatio-temporal dynamics. Previ- 
ous analysis of Refs. [14123] are extended considering two 
transverse spatial dimensions. 
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2 Equations Here E is the slowly- varying amplitude of the electric field 

within a scalar description, d cos S, d sin S are the diffusion 

Off-axis feedback loops, in which the propagation of light i i-rr i- m ■ ^ / i n\ ^ ■ ^ ^ 

^ -f o a ditrraction coethcients [d > Ojand /i is a control pa- 
is guided, are used with nonlinear phase/amplitude mod- i_ • i i i • i i 

' rameter. Devices where only the polarization evolves much 

ulators in which material evolves much more slowly than t- , i ^ • i i ■ ^ i i n i i 

faster than the electric neid are instead modelled by 



the electric field and its polarization. First examples were 
liquid crystal light valves with feedback signal propagat- 
ing through fibers whose position was easily controlled 
|9ll0lllj . These devices -in presence of a shifted feedback- where N represent the internal dynamics of the material 



dtE = gi{\E\^,N; ^i)E + de'^WlE + tc'^Eax, (3) 
dtN ^g2{\E\\N,^lN-^i), 



are modelled by equations of the type [8124) coupled with light. We consider feedback loops in which 

the time delay At is very small compared to the time scale 
of the slowly varying envelope of the field. The feedback 

(p(a;, y, t) is a real field that represents the state of the ma- can then be characterized by an amphtude < r < 1 and 

terial - for liquid crystal based devices it is related to the a phase shift = w^^*, where lol is the carrier frequency 

alignment of the molecules - at a point {x, y) and time i, and couples the field E in [x, y) with the field Eax in a 

while LpAx is evaluated at point x + Ax and time t, re- shifted point. /, gi and g2 are nonlinear complex functions 

spectively scaled with diffusion length and diffusion time, that can be derived with respect to E and N. The trivial 

The control parameter ^ is independent on x for the sake solutions of Eq. ^ and Eqs.® are E = and E = 

of simplicity, /i, /2 are real functions that can be derived N ^ N^, respectively. Interestingly, our analysis applies 

with respect to tp. In the limit of infinitely extended sys- also to equations with the more general feedback term 

tems, the homogeneous states are solutions of /i + /2 = [/i(|£'P)£^]zix, with re'*^ = /i(0). 

and their domains of existence depend upon ^ but not We are interested in the dynamics of perturbations 

upon Ax. SE ^ cxp {ujt + ik • x) around the uniform states in the 

A similar feedback is used also with active materials linear regime. The dynamics of perturbations for Eq.@ 

inside optical cavities. When the dynamics of the differ- and Eqs.® is contained in the dispersion relation 



ence of the populations of the energy levels coupled to the 
light and of the polarization of the materials is much faster 
than the dynamics of the electric field, these systems are 
described by a single equation of the type jl4l23] 



uj = j3- e'\kl + kl) + re'(0+fex4-)^ (4) 



where space is rescaled in units of ^/d. Here j3 ~ f{0) for 
Eq.© and ^ = gi(0,iVo) for Eqs.®. For class B models, 
perturbations 6N are always damped and decoupled from 
dtE ^ f{\E\^;fi)E + de'^VlE + re"f'EAx- (2) 5^; and can be ignored. 
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We note that Eq.© with (5 = 0, /3 = d^fi, r ^ \d^f2\ 
and (f) — or TT depending on the sign of |9i^/2|, is the dis- 
persion relation for the perturbations of the uniform states 
of Eq.(IT]). This shows that despite the different physical 
meaning of the variables tp and E, as well as the signifi- 
cant differences in the characteristic time scales and in the 
light-matter coupling behind Eq.(IT]) and Eq.Q or Eqs.Q, 
the dynamics of the perturbations of Eq.((T]) is a special 
case of the dynamics of the perturbations of the other two 
cases. 

From Eq.Q we find that there are unstable band of 
k = {kx, ky) with the most unstable ones given by 

VfeCj/j ~ ~{2kx cos S + rAx sin {k^Ax + 0), 2ky cos 5) 

= (0,0) (5) 
dl^ujR = -{2cos5 + rAx^ cos{kxAx + (j))) < 0, (6) 

(7) 



dl2Luii = -2cos(5 < 



The subscripts R and / refer to real and imaginary part, 
respectively. The conditions ([5][7]) ensure that the solution 
of Eq.([5l) corresponds to the perturbation with the largest 
amplification; note that Eq.© is the standard stability 
condition for diffusive equations. 

An important feature of the instability threshold is 
that it is a function of four relevant parameters, namely 
4)^5,rAx'^ , and PrAx'^, and is independent on /?/; there- 
fore increasing the shift size Ax produces on the device the 
same effect of larger gain /Sr and feedback r. A specific ef- 
fect of the nonlocality is that the relative strength of diffu- 
sion and diffraction, S, also becomes an effective parameter 
to control the threshold position. For active materials, the 
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Fig. 2. Instability threshold as a function of the shift. The con- 
trol parameter fi = Pr + 1 is represented for easy comparison 
with results in the following section. The effect of the feedback 
phase on the instability is shown: thresholds for (f) — (hor- 
izontal dark line) and increasing values — 7r/4, 7r/2, 37r/4, tt 
(lighter colors), S — 0.497r, r = 0.5. The dot corresponds to the 
parameter choice in Fig. [S] 

lowest gain and feedback thresholds are generally found in 
the purely diffractive limit {S ^ 7r/2). The threshold value 
for the scaled feedback strength rAx^ is independent on 
the sign of the refractive index (sign of 6) and increases 
with diffusion. Both /3fl and r can be increased to cross 
the laser threshold and -similarly to the case of perfect 
alignment- if the feedback is out of phase then stronger 
gain is required, as shown in Fig. [2j For vanishing shift 
(Ax = 0) this phase acts as a detuning and increases 
the threshold while for larger shift values the threshold 
tends asymptotically to the value for in-phase feedback 
{(j) ~ 0). Quite surprisingly, a misalignment lowers the 
threshold because the most unstable mode has k^j ^ 
and, in this case, the nonlocal coupling can reduce the 
feedback dephasing. Consistently with this interpretation, 
when the feedback perfectly in phase with the intracavity 
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field {(f) = 0) the most unstable mode is homogeneous with 
fix, I — and the threshold is independent on the lateral 
shift Ax. 

The dispersion Eq. ^ has in general a not vanishing 
imaginary part corresponding to a non null phase velocity 



Vp = -w/(k)' 



k|2 



(8) 



- {kl + kl)sm6 + rsm{k^ Ax + (/))]■ 



k|2' 



For the most unstable k, kc, the phase velocity is 



N /I • r Pi 2cos5 
Vp(ke) = (fc,,sm(5~|- + ^^,0), (9) 



and group velocity is 



Vg(kc) = VfcW/(kc 



{2kxc sin 5 — rZ\xcos {kxc^x + 0), 0) = 
(2fc..sinJT/^^^^^,0). (10) 



Note that both Vp(kc) and Vg(kc) are null in the direction 
orthogonal of the shift and coincide with the one dimen- 
sional values. In Eq. ([TU)) the — sign solution always sat- 
isfies Eq.([6l) and therefore corresponds to a wave-packet 
with frequency spectrum centered on a local maximum of 
the amplification. On the contrary, the -I- sign may not 
satisfy Eq.®. We remind that Eqs. pfTU]) with (5 = 0, 
(j) = Q or (j) — t: give the velocities for the diffusive systems 
of Eq. ([IJ . An example of the large variability of phase and 
group velocities at the critical wave-number are shown in 
Fig.H 

The analysis above shows that the phase velocity at 
kc and the group velocity are always parallel to the shift, 
while their sign varies. There are manifolds in the control 





Fig. 3. Variation of phase (upper plot) and group (lower plot) 
velocies, Vp,x{kc) and Vg,x{kc) respectively, depending on S and 
rAx^. Other parameters: /3 = 0.2 — i0.02 and (l> = Tr/4. 

parameter space that separate regions in which the group 
and the phase velocity of the most unstable perturbation 
have the same sign from region in which these velocities 
have opposite sign. Moreover, the real part of the disper- 
sion relation may have more than one maximum, so that a 
single perturbation may split into two wave-packets. The 
independent tunability of phase and group velocity is a 
specific feature of optical systems with two-point nonlocal- 
ity: as a matter of fact, gradient terms fix both velocities 
in the same direction of the drift, while two-point nonlo- 
cality in diffusive systems gives always opposite velocities 
at the critical wavenumber (up(kc) = — fg(kc)). 
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In the purely diffractive limit 5 Tr/2, both velocities 
are actually odd functions of kc] this symmetry is reduced 
by the effect of diffusion {\5\ < 7r/2). Therefore, even if 
for (f) — t: both +kc and —kc are unstable, from the linear 
analysis we do not expect intensity stripes above threshold 
in such optical systems. As a matter of fact, the instability 
of these two states (different traveling waves) is rather 
peculiar and opens the possibility of bistability instead of 
stripe pattern formation. 

Note that the tunability of transverse phase and group 
velocities is a general property that is valid also in the 
case in which At is of the order of the time scale of the 
slowly varying amplitude [23 . This tunability is therefore 
a rather robust and distinctive feature of two-point nonlo- 
cality with respect to models where velocities are induced 
by drift terms [53] . We remark that the parameters (p and 
6 are essential for this tunability: for this reason in the dif- 
fusive systems of Eq.(IT|) the tunability is absent and the 
phase and group velocities have always opposite sign. 

We now consider whether the perturbations with non 
null group velocities grow fast enough to occupy the entire 
system (absolute instability) or if their group velocity is 
such that the perturbations, although growing, move away 
(convective instability) . 

The analysis is performed in an infinite system eval- 
uating the asymptotic behavior of perturbations both in 
a traveling and in a fixed reference frame. Note that the 
distinction between absolutely and convectively unstable 
regimes is generally given in infinite system while the in- 
clusion of boundary effect in finite systems can change, in 
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Fig. 4. Regimes in which the homogeneous solution is stable, 
convectively unstable and absolutely unstable as a function of 
the shift and for S = 0.457r, r = 0.56, </> = 0.257r. 

some case even drastically, the instability scenario. The 
simpler example is the case of periodic boundary condi- 
tions: in this case no convective instability does exist at 
all. Another interesting related question for finite systems 
is the phenomenon of transient growth of perturbations 
observed in cases where the linear stability operator is 
non normal, i.e. does not commute with its adjoint |25) . 
The connection between convective instability and tran- 
sient growth of perturbations is still an open question |26j . 
What has been characterized is the macroscopic amplifi- 
cation of quantum noise in optical systems in this regime 

We determine absolute thresholds by evaluating asymp- 
totically the Green function: we extend analytically the 
dispersion relation to complex wavevectors k (for one trans- 
verse dimension) and find the appropriate integration paths 
in the plane k. A detailed analysis is given in [23], here 
we remark only that the for the purely diffusive systems 
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of Eq.([T|) the asymptotic evaluation of the Green func- 3 Nonlinear and stochastic spatio-temporal 



dynamics 

The general analysis of the previous section encompasses 
a broad class of non-diffractive as well as laser models. 



tion is done by closing the integration contour using only 
equiphase lines from saddle points. Furthermore, it turns 
out that the absolute threshold is determined only by 

the saddle closest to the imaginary axis. For Eq.(I2|) and 

„ . , 1 1 -,1 , ,1 , .1 Main features of the spatio-temporal dynamics of specific 

Fqs.(13|) mstead, we close with steepest descent paths a 

r, ., , ri I 1 r,i • • ■ , • systems can be anticipated from the linear stability anal- 

nnite segment [K/„,fc/j(/J of the imaginary axis contam- 

, /n , \ n 1 -,11 I ysis and we will see some examples in the case of a class 

mg all the kj with wuXU^kij > U and with k/„i, k/mj 

such that wjilO, fc/) < for fc/ < fc/^ and fc/ > fc/A/ . For 



A [Ml laser: 



these two cases, the correct determination of the absolute — —E{1 + i9 — N) + e^^V^ E + re^'^E{x + Ax y) 

thresholds requires to identify the integration paths and 
to evaluate the contribution to the Green functions of the 
saddle points that are part of it, excluding all the others. 



^ (11) 

1 + LBP ^ ' 



Summarizing, our analysis allows to distinguish three 
regimes both in diffusive {81 and laser [TT systems with 
a two-point nonlocality, as shown in Fig. |4] for a partic- 
ular parameters choice. The homogeneous (vanishing or 
non-lasing) state becomes unstable above a first threshold 
(convectively unstable regime) corresponding to positive 
dispersion Eq.Q. The absolutely unstable regime is found 
only after evaluation of the asymptotical growth of per- 
turbations, involving an integral whose approximate value 
is found with a non-trivial application of the saddle-point 

technique. This calculation is fully described in |14) and The regime of convective instability allows the amplifica- 



with ^{x, y, t) complex Gaussian white noise. Numerical 
simulations of this nonlinear model are performed with 
a second-order in time Runge-Kutta method and using 
the random number generator of Ref. [30) . The connec- 
tion with the analysis of the previous section is given by 
/3 = fj, — 1 — 10. Numerical simulations confirm the pre- 
dicted instability diagram; the wavenumbers dynamically 
selected and the velocities are well approximated by those 
obtained from the analytical analysis of the linear disper- 



sion. 



3.1 Signal splitting and interactions 



P5] and here we only stress that not monotonic threshold tion of localized light signals below the lasing threshold, 
dependence on the shift can be found (Fig. 2]) and have The direction of propagation of these signals depends on 



been also checked by numerical simulations of dynamical the group velocity Eq. (jlOp and steering can be obtained 
equations. also non-mechanically by varying the feedback loop phase. 
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Fig. 5. Spatio-temporal intensity dynamics showing the inter- 
action of signals in a ID class A laser with parameters Ax = 1, 
^ = 1.3, e = -1, S = 0.497r, r = 0.5 , = TT, e = (see Fig. 
[5]). Evolution starting from two small Gaussian perturbations 
as can be recognized in the picture at time t = 0. 

as shown in Figs. 3 and 4 of Ref. [M]. Particularly in- 
teresting is the case in which the feedback is negative, 
(f) = TT, first discussed in Rcf. T^. The recovered symme- 
try in the dispersion relations gives rise to instability of 
both positive and negative wave-vectors and, as a pecu- 
liar consequence of the two-point nonlocality, these waves 
propagate in opposite directions. This allows the system 
to operate as a signal splitter in which an initial pertur- 
bation, such as a localized spot of light, is divided into 
two counter-propagating copies. This phenomenon is ro- 
bust also in presence of noise [23] and deviations from the 
symmetry, i.e. differences between the left and the right 
propagating signals, are due to nonlinear effects |27I28| . 

Once excited, these localized signals are amplified and 
propagate in the system (even in opposite directions) with- 
out giving rise to any intensity modulation, even if the field 
is actually spatially oscillating. We will now focus on a dif- 
ferent aspect of signal control, namely signal interaction 
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when they cross each-other. As localized perturbations are 
split into counter-propagating ones it is actually possible 
to have signals crossing during temporal evolution. This 
case is represented in Fig. [5] In this example the system is 
excited in two separated points at an initial time and the 
dynamics is considered in the case of only one transverse 
dimension. We see that two of the four generated signals 
cross and locally interfere generating a transient intensity 
stripe in the center of Fig. [S] Moreover, looking on the 
left (right) side, we see that due to the different velocities 
of perturbations fronts, the trailing edge of the left signal 
is reached by the leading edge of the following signal and 
this also generates intensity modulation by interference. 
Finally, after a transient time (not shown in Fig. [5]) and 
as expected, due to the convective character of the insta- 
bility, the system evolves locally back to the homogeneous 
vanishing state. 

3.2 Patterns in one and two dimensions 

In previous works we described the spatio-temporal dy- 
namics of class A lasers assuming a one dimensional trans- 
verse geometry [14123) . Good agreement with the theoret- 
ical predictions (threshold position, pattern wave-length, 
phase and group velocity) was found when considering 
rather large systems where boundary conditions effects 
are negligible. The importance of the boundaries is par- 
ticularly evident in the convectively unstable regime. In 
Fig. [5] we show the evolution of a noise sustained structure 
for 6 — 0.457r, r = 0.56, cj) — 0.257r and super-Gaussian 
pump with maximum value /i — 0.6 in the central re- 
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Fig. 7. Temporal evolution of the intensity (upper line) and 
of the real part (lower line) of the field starting from a noisy 
50 1 00 initial condition. Parameters: Ax = 2 , = , 5 = 0.457r, r = 

0.56 , (j) = tt/A, e = 10~^. Super-Gaussian profile of /j, with 
Fig. 6. (a) Spatio-temporal evolution of the intensity field maximum value 0.6 inside the circular area. For this value of 

fi the vanishing state is convectively unstable. 





X 



for one transverse dimension and starting from a noisy initial 
condition, (b) Intensity \E{x,t = 100)^ (black line), real part 

(green dotted line) of the field, and pump profile /i (orange pattern. Oscillations appear in the field profile (phase pat- 



line). Parameters: Ax = 2, 9 = Q, 5 = 0.457r, r = 0.56 , 
(j) = -K / 4,, e = 10~^. Super-Gaussian profile of ^ with maximum 
value 0.6. 



tern given by the green line in Fig. [6}d). 

In order to give a representation of the aspect of spa- 
tial structures as they vi^ould appear in experiments with 
broad area lasers, it is interesting to consider the evolution 
gion (see Fig. From Fig. 0] and remembering that of two-dimensional field E{x,y,t). In Fig [7] we represent 
— fi — 1, we see that for these parameters the nonlas- both fields profiles and intensities for the same parame- 
ing state is convectively unstable and the field in Fig. |6^) ters as in Fig. [6l showing the aspect of a noise sustained 
has the typical incoherent profile discussed in 23 . Here structure in the convective regime. Three snapshots give 
we show that a large system is needed in order to observe an example of the dynamic character of this structure: 
noise sustained structures intensities in the system: the numerical simulations in presence of noise show an inco- 
intensity represented in Fig. ^jp is significantly high only herent traveling structure whose aspect is continuously 
in half of the system, far from the right edge of the pump changing. The importance of the boundary and system 
profile /i. In other words, the intensity growth -in this size in this case is evident as it leads to a very low inten- 
case against the shift direction- is rather slow and large sity noise sustained structure. Due to the reduced size of 
systems are needed to observe an intense noise sustained the pumped area with respect to previous ID numerical 
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Fig. 8. Section plot of the 2D intensity (black line) and real 
part (green dotted line) of the field for t — 100. The corre- 
sponding values are multiplied by a factor 10 to be compared 
with the pump profile jj, (orange line). Same parameters of Fig. 

m 

simulations (compare fj, super- Gaussians in Figs. |6] and 
[5]) the intensity profile is dominated by a smooth front 
emerging on the right side and remains very small. 

When the pump is increased this front becomes steeper 
and and a uniform intensity state occupies a broad part 
of the system. Above the absolutely unstable threshold a 
coherent phase pattern arises as shown in Fig. ^ and the 
system is in the lasing state. Starting from a noisy ini- 
tial condition, perturbations are washed away from the 
system and the intensity profile becomes intense and uni- 
form. Also the coherence of the underlying phase pattern 
increases (in this case the pattern aspect ratio is very small 
and only one oscillation is actually seen in Fig. Note 
also that the intensity profile is deformed and displaced 
against the direction of the shift, i.e. it is not exactly cen- 
tered in the pumped region (white circle in Fig. [9l panels 
at t 75). 



(-5 t-30 i~7S 

m €) 9) 



Fig. 9. Temporal evolution of the intensity (upper line) and 
of the real part (lower line) of the field for vanishing state 
absolutely unstable. Parameters: Ax = 2 , 6 — , 5 — 0.457r, 
r — 0.56 , (f) = tt/4, e — 10~^. Super-Gaussian profile of fj. with 
maximum value 1. 

The most pecuhar patterns in lasers with nonlocal 
feedback occur for (f> — tt, when the system operates as 
a signal splitter. The change in the pattern aspect when 
increasing the pump is shown in Fig. 1101 In the convec- 
tively unstable regime (two panels with /i = 1.4) there are 
now two regions of larger intensity, as perturbations with 
opposite wave-numbers travel apart, as expected. Intense 
light spots correspond to a certain coherence in the under- 
lying phase pattern, while the intensity drops down where 
there are defects in the phase stripes. Intensity reaches 
larger values on one side (here the left one, against the 
positive shift direction. Ax = 2) [14123128) . 

Noise sustained patterns for negative feedback are then 
characterized by off-axis spots and have a vanishing inten- 
sity in the central area, while crossing the absolute insta- 
bility threshold an intense and uniform profile is reached 
-after a transient- in the whole pumped region. Note that 
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Fig. 10. Spatial distributions of the intensity (left panels) 
and field real part (right panels) for numerical simulations 
with super-Gaussian pump profiles with maximum intensity 
/J, — 0.7, 1.4, and 2.1 inside the circular panels. In the exter- 
nal region the gain vanishes. Other parameters are Ax = 1, 
Q = 0.2, & = 0.497r, r = 0.5 , = TT, e = 10"^. Simulations are 
started from noise initial condition. 

the phase pattern has larger wave-numbers in the noise 
sustained structure. After a longer transient (not shown 
in the picture) the stripe in the absolutely unstable regime 
(Fig. [To] for fi — 2.1) becomes orthogonal to the shift di- 
rection, as theoretically predicted. 



Two-point nonlocality can be easily induced in optical de- 
vices through a displaced feedback loop and leads to dis- 
tinctive effects opening interesting possibilities in experi- 
ments in extended devices. We have presented its analyti- 
cal characterization through linear approximations, allow- 
ing to predict instability convective and absolute thresh- 
olds, unstable wave-numbers and velocity of drifting pack- 
ets, in a broad class of two-dimensional nonlinear systems. 
Many differences with respect to the most studied drift, 
modelled by a gradient term, make this two-point nonlo- 
cality an interesting and versatile tool in view of experi- 
ments. In particular, such nonlocality opens significantly 
larger windows of control parameters where the system 
output state is sustained by the presence of noise instead 
of the dynamics (convective regime). In such regime the 
system can be locally excited generating a signal that is 
amplified during propagation. Note that the amplification 
here considered occurs once the (initial) local perturba- 
tion is removed, at difference of laser amplification above 
transparency, generally considered in presence of a contin- 
uous signal injection. 

Another distinctive feature of displaced feedback and 
two-point nonlocality is the possibility to independently 
tune phase and group velocity, while gradient terms fix 
both velocities in the same direction of the drift. In optical 
systems such velocities have been shown to be either par- 
allel or opposite and can be tuned even non-mechanically 
through the feedback phase ip (without touching the ex- 
periment alignment). Such large tunability is a charac- 
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teristic of optical systems, while two-point nonlocality in lytical predictions for infinite systems are in quantitative 

diffusive systems gives always opposite velocities at the agreement with simulations for relatively large systems, 

critical wavenumber (wp(kc) — — Ug(kc)). Another distinc- while a finite system effects gives important deviations in 

tive feature of two-point nonlocality in optical systems is small systems not only when few oscillations of the phase 

the possibility to control the feedback phase to operate the patterns are present but also in the cases in which smooth 

laser as a signal splitter (for (j) — t:). Then an initial per- fronts modulate the light intensity profile, 
turbation is split in two copies ampfifying (even when the Funding from FISICOS (FIS2007-60327) and CoQuSys 

injection is removed) during their propagation in opposite (200450E566) is acknowledged, 
directions through the broad area device. 
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